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ABSTRACT 

We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by 
comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use 
the c2d cata l ogs to separate 108 s tarles s fro m 92 protoste l lar co res in the 1.1 mm core samples from 
lEnoch et al.1 (|2006l ). lYoung et al.l (|2006ft . and lEnoch et~afl (|2007f >. A comparison of these populations 
reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but 
larger sizes and lower densities on average than protostellar cores, with sizes that suggest density 
profiles substantially flatter than p oc r~ 2 . By contrast, starless cores in Serpens are compact and 
have lower masses than protostellar cores; future star formation will likely result in lower mass objects 
than the currently formin g protostars. Compariso n to dynamical masses estimated from the NH3 
survey of Perseus cores by iRosolowskv et all (j2008f ) suggests that most of the starless cores are likely 
to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution 
includes 108 cores and has a slope of a = —2.3 ± 0.4 for M > 0.8 M Q . This slope is consistent with 
recent measurements of the stellar initial mass function, providing further evidence that stellar masses 
are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency 
of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud, thus 
the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5 x 10 5 
years, with a total uncertainty of a factor of two. Such a short lifetime suggests a dynamic, rather than 
quasi-static, core evolution scenario, at least at the relatively high mean densities (n > 2 x 10 4 cm~ 3 ) 
to which we are sensitive. 

Subject headings: stars: formation — ISM: clouds — ISM: individual (Perseus, Serpens, Ophiuchus) 
- submillimeter - infrared: ISM 



1. INTRODUCTION 

Dense prestellar cores, from which a new generation 
of stars will form, represent a very early stage of the low 
mass star formation process, before collapse results in the 
formation of a central protostar. The mass and spatial 
distributions of these prestellar cores retain imprints of 
their formation process, and their lifetime is extremely 
sensitive to the dominant physics controlling their forma- 
tion. It is well established that most of the star formation 
in our Galaxy occurs in clusters and groups within large 
molecular clouds (e.g., lLada fc Ladal I2003L and refer- 
ences therein). Molecular clouds are k nown to be turbu- 
lent, with supersonic line-widths (e.g. lMcKee &; Zweibell 
1992), and to have complex magneti c fields that are likely 
important to the cloud physics (e.g.. lCrutcherll 999). Un- 
derstanding the properties of prestellar cores on molec- 
ular cloud scales, and how they vary with environment, 
provides insight into the global physical processes con- 
trolling star formation in molecular clouds. 
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For example, in the clas sic paradigm of magnetically 
dominated star formation (|Shu et al.lll987l ). the collapse 
of cores occurs very slowly via ambipolar diffusion and 
starless cores sho uld be long liv ed, with lifetimes of or- 
der £ad ~ lOtff (Nakano 1998), where iff is the free- 
fall timescale (tg ~ 10 5 yr for n ~ 10 5 cm -3 , where 
tff depends on the mean core density n: tg cx n~ - 5 ). 
Alternatively, if molecular cloud evolution is driven pri- 
marily by turbulence, over-dense cores should collapse 
quic kly, on approximately a dynamical timescale, 1 — 
2iff (iBallesteros-Paredes et al.ll2003llMac Low k, Klessenl 
|2004| ). Thus, the lifetime of prestellar cores should be a 
strong discriminator of core formation mechanisms. Pub- 
lished measurements of the prestellar core lifetime vary 
by two orders of magnitud e, however, from a few times 
10 5 years to 10 7 years (see Ward-T hompson et aTll2007l 
and references therein). 

While the properties of star-forming cores depend 
strongly on the physical processes leading to their for- 
mation, core initial conditions in turn help to determine 
the evolution of newly formed protostars. One of the 
most important diagnostics of initial conditions is the 
mass distribution of prestellar cores. In addition to be- 
ing a testable prediction of core formation models, a 
comparison of the core mass distribution (CMD) to the 
stellar initial mass function (IMF) may reveal what pro- 
cess is responsibl e for determining stellar masses (e.g., 
iMever et al.ll2000D . 

A number of recent studies have found observa- 
tional evidence that the shape of the IMF is directly 
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tied to the core fragmentation process (iTesti fc Sargent 
ll998tlMotte et al.lll998HOnishi et al1l2002l ). lAlves et al. 
1)20071 ) find a turnover in the CMD of the Pipe Neb- 
ula at M ~ 3x th e turnover in the Trapezium IMF 
(jMuench et al.ll2002t ). and suggest that the stellar IMF 
is a direct product of the CMD, with a uniform core- to- 
star efficiency of 30% ± 10%. The mean particle densities 
of the extinction-identified cores in Pipe Nebula study 
(5 x 10 3 — 2 x 10 4 cm~ 3 ) are lower than those of typi- 
cal cores traced by dust emission (2 x 10 4 — 10 6 cm~ 37 ; 
lEnoch et al.ll2007ti , however , and the cores may not be 
truly prestellar. In Orion, iNutter fe Ward-Thompson] 
(|2007Tl find a turnover in the CMD of starless SCUBA 
850/im cores at ~ 1.3M Q . T hose authors re late this 
turnover to a down-turn in the iKroupal (|2002t ) IMF at 
~ O.IMq, and infer a much lower core-to-star efficiency 
of 6%. 

Large samples of prestellar cores are important for fur- 
ther addressing these problems, as is a more reliable sep- 
aration of prestell ar, protostellar, and unbo und starless 
cores. We follow iDi Francesco et ail (|2007T > in defining 
starless cores as low mass dense cores without a compact 
internal luminosity source, and prestellar cores, at least 
conceptually, as starless cores that are gravitationally 
bound and will form stars in the future. For millimeter 
cores containing a compact luminous internal source (i.e., 
an em bedded protostar) we follow IDi Francesco et ail 
(2007) in terming these protostellar cores, regardless of 
whether the final object will be stellar or sub-stellar in 
nature. Unlike protostellar cores, which are internally 
heated by the embedded source as well as externally by 
the interstellar radiation field (ISRF), starless cores are 
heated only externally by the ISRF, w ith decreasing tem - 
peratures toward the core center (e.g.. lEvans et al.ll2001f ). 

Molecular line or extinction surveys often trace rela- 
tively low density material (10 3 — 10 4 cm -3 ), leaving the 
possibility that such cores may never collapse to form 
stars. In addition, most previous studies base the iden- 
tification of protostellar versus starless cores on near- 
infrared data, which is not sensitive to the most embed- 
ded protostars, or on low resolution and poor sensitivity 
IRAS maps. The first issue can be remedied by using 
millimeter or submillimeter surv eys; (sub)mm emission 
trace s dense (n > 2 x 10 4 cm~ 3 : Ward-Thompson et"ai1 
11994 lEnoch et all l2007t ) material, and detection at 
(sub)mm wavelengths tends to correlate with other in- 
dicati ons of a prestellar natur e, such as inward mo- 
tions (|Gregersen fc Evansl 12000) . Spitzer provides sig- 
nificant progress on the second issue, with substantially 
superior reso l ution and sensitivity (0.01L© at 260 pc; 
lHarvev et al.l l2007af ) compared to IRAS, making the 
identification of prestellar cores much more secure. Re- 
cent exam ples of identifyin g starl ess cor es include the 
studie s of Dorgensen et al.l (|2007h and lHatchell et al.1 
l|2007h . which utilize SCUBA 850 /xm surveys and Spitzer 
data to distinguish starless from protostellar cores in the 
Perseus molecular cloud. 

The combination of (sub)mm studies with molecular 
line observations of dense gas tracers, which yield a gas 

7 The mean particle density of Bolocam cores is given by n = 
3M/(4irR(im[i), where M is the core mass, R the radius, fi = 2.33 
is the mean molecular weight per particle, and mu the mass of 
Hydrogen. 



temperature and line- width, is a powerful method for de- 
termining the mechanical balance of starless cores. While 
our preliminary operational definition of a prestellar core 
will be a starless core that is detected at submillime- 
ter or millimeter wavelengths, we will examine this is- 
sue more closely in § |4] Comparison of our data to 
molecular line observations, such as the recent GBT NH3 
(2,2) and (1,1) survey o f Bolocam cores in Perseus by 
iRosolowskv et al.l (|2008|) provides a more robust method 
of estimating whether cores are gravitationally bound. 

We have recently completed large continuum surveys 
at A = 1.1 mm of the Perseus, Ophiuchus, and Ser- 
pens molecular clouds using Bolocam at the Caltech 
Submillimeter Observatory (CSO). We mapped 7.5 deg 2 
(140 pc 2 at our adopted cloud distance of d = 250 pc) 
in Perseus, 10.8 deg 2 (50 pc 2 at d = 125 pc) in Ophi- 
uchus, and 1.5 deg 2 (30 pc 2 at d = 260 pc) in Serpens 
with a resolution o f 31" (|Enoch et aT]|2006t lYoung et all 
I2006t lEnoch et al.ll2007l hereafter Papers I, II, and III, 
respectively). Millimeter emission traces the proper- 
ties of starless cores and protostellar envelopes, includ- 
ing core sizes, shapes, masses, densities, and spatial 
distribution. The 1.1 mm Bolocam surveys are com- 
plemented by large Spitzer Space Telescope IRAC and 
MIPS maps of the same clouds from the "From Molecu- 
lar Cores to Planet-forming D isks" Spitzer Lega cy pro- 
gram ("Cores to Disks" or c2d: lEvans et al.ll2003l ). Com- 
bining these data sets with the 2MASS survey provides 
wavelength coverage from 1.25 — 160/im, and enables us 
to reliably differentiate starless cores from those that 
have already formed embedded protostars. The Spitzer 
c2d IRAC and MIPS surveys of each cloud are de- 
scribe d in detail in f jorgenscu et al. (2006), Harve v~et al.1 



(2006), R ebull et al 
lHarvev et al.l (12007, 



2007h . iPadgett et al.l (|2007h . and 



In Paper III, we looked at how the global molecu- 
lar cloud environment influences the properties of star- 
forming cores, by comparing the 1.1 mm core populations 
in Perseus, Serpens, and Ophiuchus. In a companion pa- 
per to this work (Enoch et al. 2008, in prep), we use 
the comparison of 1.1 mm and Spitzer data to study the 
properties of embedded protostars. In particular, we ex- 
amine the bolometric temperatures and luminosities of 
Class and Class I protostars in Perseus, Serpens, and 
Ophiuchus, the lifetime of the Class phase, and accre- 
tion rates and history for the early protostellar phases. 
Here we use the combination of Bolocam 1.1 mm and 
Spitzer c2d surveys to probe the initial conditions of star 
formation on molecular cloud scales, and how the prop- 
erties of starless cores differ from cores that have already 
formed protostars. 

In we describe the identification of protostellar 
cores, including the combination of 1.1 mm and Spitzer 
infrared (IR) data (§ 12. ip . identification of candidate pro- 
tostars based on their mid- and far-infrared properties 
( ^2.2p . and the basis on which we determine association 
between 1.1 mm cores and candidate protostars 12. 3p . 
The resulting starless and protostellar 1.1 mm core pop- 
ulations for each cloud are compared in including core 
sizes and shapes ( £13. ip . masses and densities ( £|3.2p , dis- 
tribution in mass versus size ( §3.3|) . relationship to cloud 
column density ( ^3 .4|) . and spatial clustering (<JT5|). We 
calculate the dynamical mass of starless cores in Perseus 
using NH3 observations to determine if the cores are truly 
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prestellar (§3]). We combine the three clouds to produce 
the prestellar core mass distribution (CMD), which is 
discussed in relation to the stellar initial mass function 
in <j5] Finally, in § [6] we estimate the lifetime of the 
dense prestellar core phase and discuss implications for 
star formation theory. 

2. SEPARATING STARLESS AND PROTOSTELLAR CORES 

To study the initial conditions of star formation as 
traced by prestellar cores, we first must differentiate cores 
without an internal source of luminosity (starless cores) 
from those with an embedded self-luminous source (pro- 
tostellar cores). Protostellar cores will have lost some 
mass due to accretion onto the embedded protostar, and 
may be otherwise altered, so that they are no longer rep- 
resentative of core initial conditions. Starless and proto- 
stellar cores can be differentiated using the Spitzer c2d 
surveys, by identifying infrared sources that may be as- 
sociated with a given core. Such candidate protostars 
are typically visible as point-like objects in the near- to 
mid-infrared data. In the following sections we describe 
the merging of the millimeter and infrared data and the 
criteria used to determine which cores are protostellar. 

2.1. Combining Bolocam and Spitzer c2d Data 

Spitzer IRAC and MIPS maps from the c2d Legacy 
program cover nearly the same area as our Bolocam 
1.1 mm maps of Perseus, Serpens, and Ophiuchus. Both 
Bolocam and Spitzer maps were designed to cover down 
to a visual extinction of Ay > 2 mag in Perseus, 
Ay > 3 mag in Ophiu chus, and Ay > 6 mag in Ser- 
pens ([Evans et al.| [2003y The actual overlap in area be- 
tween Bolocam and IRAC maps is shown in Figure 1 
of Papers I, II and III for Perseus, Ophiuchus, and Ser- 
pens, respectively. Catalogs listing c2d Spitzer fluxes of 
all detected sources in each of the three clouds, as well 
as near-infrared fluxes for sources that also appear in 
the 2MA SS catalogs, are av ailable through the Spitzer 
database (|Evans et al.l l2007). Thus, we have wavelength 
coverage from A = 1.25 to 1100 /zm, utilizing 2MASS 
(A = 1.25, 1.65, 2.17 fim), IRAC (A = 3.6, 4.5, 5.8, 
8.0 fim), MIPS (A = 24, 70, 160 pm), and Bolocam 
(A = 1.1 mm) data. Note that 160 /jm flux measure- 
ments are not included in the c2d delivery catalogs due 
to substantial uncertainties and incompleteness, but arc 
utilized here and in Enoch et al. (2008, in prep.) when 
possible. 

Basic data papers describe the processing and anal- 
ysis of the Spitzer IRAC and MIPS maps of Perseus, 
Serpens, and Ophiuchus, and present general proper- 
ties of the sources in each clo ud such as colo r -color 
and color- ma gnitud e diagrams (iJorgensen et all 120061: 
Harvey et all 120061: iRebull et all 120071 : lHarvev et all 
2007al) . In addition, the young stellar obje ct (YSO) pop- 
ulation in Serpens is discussed in detail bv lHarvev et alJ 
(|2007bf) . Here we are most interested in very red sources 
that are likely to be embedded in the millimeter cores 
detected with Bolocam. For the following we will use 
the term "candidate protostar" in ge neral to encompass 
candidate Class and Class I objects (jAndre et al.lll993l : 
lLada fc Wilkmel Il98l . although more evolved sources 
may be included in this sample as well. A more detailed 
study of the properties of the candidate protostars them- 
selves is carried out in a companion paper (Enoch et al., 



in prep). 

In Figure [T] we show the result of combining Spitzer 
and Bolocam data for a few cores in each cloud. Three- 
color (8, 24, 160 (im) Spitzer images are overlaid with 
1.1 mm Bolocam contours, and symbols mark the posi- 
tions of all identified 1.1 mm starless cores ("x"s) and 
protostellar cores (squares). Bolocam IDs of the central 
sources, from Table 1 of Papers I— III, are given at the 
top of each image. Red 160 fim images are often af- 
fected by saturation, pixel artifacts (bright pixels), and 
incomplete coverage. Saturation by bright sources af- 
fects many pixels in a given scan; this and incomplete 
sampling accounts for the striped appearance of the red 
images, particularly in Ophiuchus. 

Note the lack of infrared point sources near the cen- 
ter of the starless 1.1 mm cores, whereas protostellar 
1.1 mm cores are clearly associated with one or more 
Spitzer sources. Apart from these common traits, both 
starless and protostellar cores display a wide range of 
properties. They may be isolated single sources (e.g. 
Per-Bolo 62, Oph-Bolo 42), associated with filaments or 
groups (Per-Bolo 74, Ser-Bolo 20), or found near very 
bright protostars (Per-Bolo 57). Although most of the 
starless cores are extended, some are quite compact (e.g., 
Per-Bolo 57), barely resolved by the 31" Bolocam beam. 
Some starless cores are distinguished in the Spitzer bands 
by their bright 160 /im emission. Per-Bolo 62 is not de- 
tectable at 24 or 70 /im, but emerges as a diffuse source 
at 160 fim, indicating a cold, extended core that closely 
mirrors 1.1 mm contours. Other starless cores stand out 
as dark "shadows" in the shorter wavelength bands. Ser- 
Bolo 18 and Oph-Bolo 26 are two examples of such cores; 
the dark shadows against bright 8 and 24 /an emission 
are suggestive of dense cores obscuring background neb- 
ular emission. Again, contours at 1.1 mm closely trace 
the Spitzer short-wavelength shadows. 

Although the sky coverage of the IRAC, MIPS and 
Bolocam maps overlaps nearly perfectly for our purposes, 
there is a small portion of the Serpens 1.1 mm map that 
is not covered by the 70 (im map. In addition, the 160 /jm 
maps are often saturated near bright sources and in re- 
gions of bright extended emission, such as near bright 
clusters of sources. 



2.2. Identifying Candidate Protostars 

The identification of candidate protostars is discussed 
in more detail in Enoch et al. (2008, in prep); here we 
briefly describe the sample used to identify protostellar 
cores. The first criteria used to select candidate proto- 
stars from the c2d source catalogs is based on the source 
"class." All sources in the c2d database are assigned a 
class parameter based on colors, magnitudes, and stel- 
lar SEP fits, as dis cussed i n the c2d Delivery Do cument 
(jEvans et al.ll2007t ) and in lHarvev et~aT1 (|2007aD . Class 
parameters include "star" , "star+disk" , "YSOc" (young 
stellar object candidate), "red", "rising", "Gale" (galaxy 
candidate), etc. Protostar candidates will generally be a 
subset of YSOc sources, but some of the most embedded 
may also be assigned to the "red" class if they are not 
detected in all four IRAC bands. 

A preliminary sample is formed from sources classi- 
fied as "YSOc" or "red." We impose a flux limit at 
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Fig. 1. — Three color Spitzer images (8.0, 24, 160 fim) of selected starless and protostellar cores in Perseus, Serpens, and Ophiuchus. 
Bolocam 1.1 mm contours are overlaid at intervals of 2, 3, 5, ...15, 20. ..35 a, where la is the mean rms noise in each 1.1 mm Bolocam map 
(15 mjy in Perseus, 10 mjy in Serpens, and 25 mjy in Ophiuchus). Positions of starless cores are indicated by an "x", protostellar cores by 
a square, and the Bolocam ID of the centered core is given. Note the lack of infrared point sources near the starless core positions. Starless 
cores are most notable in the infrared by their extended 160^tm emission or dark shadows at shorter wavelengths, while protostellar cores 
are clearly associated with Spitzer point sources. Saturation, image artifacts, and incomplete sampling causes the striped appearance of 
some of the red 160 fim images. 

A = 24/im of S^/im > 3 mjy 8 , high enough to eliminate 
most extragalactic interlopers and sources with SEDs 

8 Note that this limit is considerably higher than the S24/xm ^ 
0.7 mjy limit imposed by Harvey ct al. (2007b), together with sev- 
eral other criteria, to help eliminate galaxies. 



that are clearly inconsistent with an embedded nature 
(e.g., sources with S^/xm < Ss/jm), but low enough to 
include most deeply embedded known protostars. In ad- 
dition, we include any 70 /im point sources that are not 
classified as galaxy candidates ("Gale"). In each cloud, 
a few known deeply embedded sources that have strong 
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Fig. 2. — Left: Distribution of the positional offset from each protostar candidate in Perseus to the nearest 1.1 mm core centroid. 
Distances are in units of the core FWHM size (8i mm ), and circles enclose sources within 0.5 and 1.0 #i mm of a core position. Pointing 
errors in the Bolocam maps are typically 10" or less. Protostar candidates include all "YSOc", "red" and 70 fim sources, while candidates 
with bolometric temperatures T^ i < 300 K (Enoch et al. 2008, in prep) are shown in bold. These cold objects, which are most likely to be 
most embedded, are primarily located within 1.0 9\ mm of a core centroid. Right: Expected distribution if the IR sources were distributed 
randomly over the Bolocam map area (i.e. if they were not spatially correlated with the cores). Our definition that a given 1.1 mm core is 
protostellar if there is a candidate protostar within 1.0 9\ m m of the core position will result in a few false associations in each cloud. 



70 /<m emission but very weak 24 /im emission (e.g., 
HH211 in Perseus) are recovered by this last criteria, as 
are a few very bright sources that are saturated at 24 \m\ 
(these are often classified as "rising"). The Spitzer c2d 
surveys are complete to young objects with luminosi- 
ties as low as 0.05L (T) (Dunham et al. 2008, in prep.; 
lHarvev et al"1l2007a[) . and we are unlikely to be missing 
any protostellar sources down to this level (Enoch et al. 
2008, in prep.; Dunham et al. 2008, in prep.). 

2.3. Determining association between cores and 
candidate protostars 

We next determine which of the 1.1 mm cores are as- 
sociated with a candidate protostar. The positional off- 
set from each candidate protostar to the nearest 1.1 mm 
core centroid position, in units of the core full-width 
at half-maximum (FWHM) size, is plotted in Figure [2] 
(le ft). Note that this a nalysis is similar to Figure 2 
of lJ0rgensen et all (|2007f ). Large circles enclose candi- 
date protostars that are located within 0.5 x 9i rnm and 
1.0 x 0\ mm of a 1.1 mm core position, where 9i mm is the 
angular FWHM size of a given core (the measurement 
of core sizes is described in § I3.1[) . Protostellar sources 
with bolometric temperature Tboi < 300 K (Enoch et al. 
2008, in prep.) are indicated by bold symbols. In gen- 
eral, the coldest objects, those expected to be embedded 
within millimeter cores, arc located within 1.0 x 9\ mm 
of a 1.1 mm core position. Thus we define protostellar 
cores to be those cores which have a candidate protostar 
located within 1.0 x B\ mm of the core center. Note that 
protostellar cores are defined by any candidate protostar 
within 1.0 x 0i mm , not only those with T^i < 300 K. 

The right panel of Figure demonstrates what we 
would expect if the same sample of Spitzer sources was 



distributed randomly over the Bolocam map area. Al- 
though there are approximately the same number of 
sources in the random sample as in the candidate pro- 
tostar sample, most do not appear on the plot because 
they are located much farther than 2 x 9i mm from the 
nearest core. There are 5 sources from the random sam- 
ple located within 1.0 x 9\ mm of a core position; thus we 
can expect a few false associations between 1.1 mm cores 
and candidate protostars based on this criteria. Indeed, 
some associations between cores and Spitzer sources are 
probably just projections on the sky. We are especially 
skeptical of associations where the 1.1 mm core flux den- 
sity is much higher than the 70 /im flux. Using 0.5 x 0\ m m 
would be a more restrictive choice, but might mis-identify 
some protostellar cores as starless. As we do not want to 
contaminate our starless core sample with more evolved 
sources, we adopt the more conservative criteria. Thus 
the number of starless cores in each cloud is likely a lower 
limit to the true value; using 0.5 x 0\ mm would result in 
3 more starless cores in Perseus, 2 in Serpens, and 4 in 
Ophiuchus. 

3. COMPARING THE STARLESS AND PROTOSTELLAR 
1.1 MM CORE POPULATIONS 

Tables |4] and [5] list the Bolocam identifications, posi- 
tions, and peak flux densities (from Papers I— III) of star- 
less and protostellar cores in each cloud. We find a total 
of 108 starless and 92 protostellar cores in the three cloud 
sample. Cores are identified based on a peak flux density 
at least 5 times the local rms noise level, and positions 
are determined using a surface brightness-weighted cen- 
troid, as described in Paper III. The Bolocam maps of 
Perseus, Serpens, and Ophiuchus are shown in Figure [31 
with the positions of starless and protostellar cores indi- 
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cated. The average la rms is 15 mJy beam -1 in Perseus, 
10 mJy beam -1 in Serpens, and 25 mJy beam -1 in Ophi- 
uchus. 

General core statistics, including the number of star- 
less (A^sl) and protostellar (iVpg) cores in each cloud, as 
well as the ratio N^h/Npg, are given in Table [TJ Note 
that the number of starless and protostellar cores are 
approximately equal in each cloud (Nsl/Nj>s = 1.2 in 
Perseus, 0.8 in Serpens, and 1.4 in Ophiuchus), a fact 
that will be important for our discussion of the starless 
core lifetime in $6] The last column of Table Q] gives the 
number of individual 1.1 mm cores that are associated 
with more than one candidate protostar (each located 
within 1 x 9\ mrn of the core position). There are 13 such 
"multiple" protostellar sources in Perseus (24% of the 
protostellar core sample), 11 in Serpens (55%), and 3 in 
Ophiuchus (17%). In general there are two or three can- 
didate protostars associated with each "multiple" core, 
with the exception of one core in Serpens (5 candidate 
protostars). Throughout this work multiple protostellar 
cores are treated as single objects. 

We now compare the physical properties of the starless 
and protostellar core populations in each cloud, with two 
primary goals. Isolating a starless sample allows us to 
probe the initial conditions of star formation, and differ- 
ences between the starless and protostellar core samples 
are indicative of how the formation of a central protostar 
alters core properties. In the following sections we follow 
the methodology of Paper III, examining the sizes and 
shapes of cores, their peak and mean densities, distribu- 
tions of core mass versus size, spatial clustering proper- 
ties, and relationship to the surrounding cloud column 
density. 

3.1. Sizes and Shapes 

Source angular FWHM sizes (9 me as) are measured by 
fitting an elliptical Gaussian after masking out nearby 
sources using a mask radius equal to half the distance 
to the nearest neighbor (see Paper II). The angular de- 
convolved core size is the geometric mean of the decon- 
volved minor and major angular FWHM sizes: 9 dec = 

y/9d,maj 9d,min, where 9d — \/ 8tneas ~ @mb an< ^ ® mb ~ 

31" is the beam FWHM. Deconvolved sizes for starless 
and protostellar cores are given in Tables [4] and El respec- 
tively, and the size distributions are plotted in Figure 01 
As discussed in Paper III, the measured size does not 
necessarily represent a physical boundary, but rather is 
a characteristic scale that depends on the linear reso- 
lution and intrinsic source density profile. For sources 
with power law density profiles, which do not have a well 
defined size, 9dec/9 m b is independent of dist ance and sim- 
ply r elated to the index of the power law (| Young et al.l 
2003). According to the correlation betw een 9dec/9 m b 
and d ensity power law exponent p found bv lYoung et al.l 
(2003!), a mean 6 dec /6 mb value of 50"/31" = 1.6 for pro- 
tostellar cores in Perseus implies an average power law 
index of p ~ 1.4 to 1.5. Many well-known protostel- 
lar sources have been found to have envelopes consis- 
tent with power law density profiles, as determined by 
high-resolution imaging combined with radiative trans- 
fer modeling. Our inferred average index for Perseus 
protostellar cores (p — 1.4 — 1.5) is consistent with the 
meanp ^1.6 from radiative transfer modeling of Class 



and Class I envelopes (jShirlev et al.l l2002t I Young et al.1 
2003). For reference, a singular isothermal sphere (SIS) 
has p = 2 (9dec/9 m b ~ 0.9), and the profile expecte d for 
a free-falling envelope i s p — 1.5 (9dec/9 m b ~ 1-6) (|Shul 
[19771 lYoung eTaIll2003h . 

Starless cores in Perseus are larger on average than pro- 
tostellar cores (Figure 13} ; the starless distribution is rel- 
atively flat, with a few barely resolved cores, and several 
larger than 3 9 m b- The Student's T-test confirms that the 
mean sizes of starless and protostellar cores are signifi- 
cantly different (significance of the T statistic is 2x 10 -5 ). 
A mean 9 de c/9mb of 2.2 for starless cores in Perseus would 
imply an extremely shallow mean power law index of 
p ~ 1.1, and the maximum value (9d e c/9 m b ~ 3.5) would 
correspond to p < 0.8. Other flat tened profi l es, such 
as th e Bonnor-Ebcrt (BE) sphere (jEbertJ fl955t iBonnorl 
Il956| ). could also produce large 9 de c/9mb values. A BE 
profile with a central density of 10 5 cm -3 and an outer 
radius of 6 x 10 4 AU would correspond to 9 d ecl^mb = 2.0 
at the distance of Perseus. There is significant observa- 
tional evidence that many starless cores do indeed look 
like BE spheres (e.g..Uohnstone et aL 2000; S hirley et al.1 
120001: lAlves et al.ll200lHEvans et al.ll2001h . We conclude 
that very large starless cores are more consistent with BE 
spheres or other flattened density profiles (p < 1) than 
with the classical pre-collapse SIS (p = 2). 9 

Despite their size differences, there is little difference in 
the axis ratios of starless and protostellar cores in Perseus 
(Figure [5]). Both populations are slightly elongated on 
average. The axis ratio is defined at the half-max con- 
tour, using deconvolved sizes: 9d,maj/Qd,min- Values for 
individual cores are given in Tables |4| and [5l Starless and 
protostellar cores have mean axis ratios of 1.7 ± 0.9 and 
1.8 ± 0.9. Standard deviations quoted here and in the 
rest of are dispersions in the sample, not errors in the 
mean. Monte Carlo tests (Paper I) indicate that cores 
with axis ratios less than 1.2 should be considered round. 

In contrast to Perseus, starless cores in Serpens are 
no larger than protostellar cores. Mean values for both 
(9 d ec/0mb ~ 1.5) correspond to an average power law 
index of p = 1.5, similar to that found for protostellar 
cores in Perseus, and to typical radiative transfer mod- 
eling results for protostellar envelopes. Starless cores in 
Serpens may be more elongated than protostellar cores 
(mean axis ratios of 1.7±0.5 and 1.4±0.3, respectively), 
but the difference is not statistically significant. 

In Ophiuchus there is no measurable difference between 
the starless and protostellar populations. The starless 
and protostellar samples have similar mean sizes, and 
both display a bimodal behavior (Figure [5). The lower 
peak is similar to the single peak seen in Serpens at 
9 dec — 1 — 2 9 m b, while the smaller upper peak is at 
sizes of 9 dec = 3 — 4 9 m b, comparable to the largest star- 
less cores in Perseus. The mean 9 de c/9mb for protostellar 
cores (1.9) corresponds to an average power law index of 
p = 1.3. Both starless and protostellar cores in Ophi- 
uchus appear fairly round, with the mean axis ratios of 
1.3 ± 0.2 and 1.4 ± 0.5, respectively (Figure 0). As dis- 

9 As the pressure-truncated boundary radius of the BE solu- 
tion has not been observationally verified, the resemblance between 
starless cores and BE profiles is necessarily limited to the region 
inside the boundary radius. Also note that a sufficiently centrally 
condensed BE sphere will be indistinguishable from a power law 
on the scales to which we are sensitive. 



Prestellar cores in Perseus, Serpens, and Ophiuchus 



7 



TABLE 1 

Statistics of 1.1 mm cores in the three clouds 



Cloud 


Ntotal 1 


N SL 2 


N PS 3 


Nsl/Np S 


N PS (mult) 4 


Perseus 


122 


67 


55 


1.2 


13 


Serpens 


35 


15 


20 


0.8 


11 


Ophiuchus 


43 


26 


17 


1.5 


3 



1 Total number of identified 1.1 mm cores. 2 Number of starless 1.1 mm cores, i.e., cores that do not have a protostar candidate located 
within 1.0 X 9i mm of the core position. 3 Number of protostellar cores. 4 Number of protostellar cores that are associated with more than 
one candidate cold protostar (each within 1.0 X 8\ mm of the core position). 



cussed in Paper III, larger axis ratios in Perseus and Ser- 
pens may be at least partly an effect of the lower linear 
resolution in those clouds compared to Ophiuchus (e.g. 
blending) . 

Were we to calculate linear sizes, cores in Ophiuchus 
would appear smaller by nearly a factor of two compared 
to Perseus and Serpens, given the smaller distance to 
Ophiuchus. This is primarily a systematic effect of the 
linear resolution, however. In Paper III we found that 
convolving the Ophiuchus map with a larger beam to 
match the linear resolution of Perseus and Serpens pro- 
duced larger linear core sizes by nearly a factor of two, 
but similar measured Odec/^mb values. Thus we focus 
here on angular sizes only. 

To summarize, protostellar cores in all three clouds 
have mean sizes consistent with power law density pro- 
files with an average index p = 1.3 — 1.5. Starless cores 
in Perseus are significantly larger on average than proto- 
stellar cores, suggestive of BE spheres or other shallow 
density profiles. Starless cores are quite compact in Ser- 
pens, while both starless and protostellar cores in Ophi- 
uchus display a bimodal distribution of sizes, with a few 
very large cores. The deficit of cores with 6dec/0mb > 2 
in Serpens and Ophiuchus as compared to Perseus may 
be related to the general lack of isolated sources in those 
clouds; the measured size of a core is limited by the dis- 
tance to the nearest neighboring source, so sizes will tend 
to be smaller in crowded regions (e.g. L 1688 Ophiuchus) 
than for isolated sources. 

3.2. Core Masses and Densities 

Core total masses, given in Tables |4] and [5j are calcu- 
lated from the total 1.1 mm flux, Sx.Xmm : 

d 2e S 

where al is the cloud distance, B u is the Planck function 
at dust temperature To, and Kx.Xmm = 0.0114 cm 2 g _1 
is the dust opacity. The total 1.1 mm flux density is 
integrated in the largest aperture (30" — 120" diameters 
in steps of 10") that is smaller than the distance to the 
nearest neighboring source. We assume that the dust 
emission at A = 1 . 1 mm is optically thin, and that To and 
/ti.imm are independent of position within a core. The 
v alue of Kx.imm is interpolated from Table 1 column 5 
of lOssenkopf fc Hennind (|1994h for dust grains with thin 
ice mantles, and includes a gas to dust mass ratio of 100. 
A recent measurement of k based on near-IR data and 
450/xm an d 850/zm SCUBA m aps yields Kx.Xmm = 0.0088 
cm e" -1 ([Shirley et al.ll2007l) . which would increase our 
masses by a factor of 1.3. 

The value of To should depend on whether a core is 
starless or has an internal source of luminosity, so we as- 
sume a slightly higher temperature for protostellar cores 



(Td = 15 K) than for starless cores (T D = 10 K). For 
dense regions without internal heating, the mean temper- 
ature is about 10 K, warmer on the outside and colder 
on the inside (|Evans et alj 12001). A recent NH3 survey 
of the Bolocam cores in Perseus confirms that the me- 
dian kinetic temperature of starless cores in Perseus is 
10.8 K (see Schnee et al. 2008, in prep, and §gj. Our 
assumed value of 15 K for protostellar cores is the average 
isothermal dust temperature found from ra diative trans- 
fer models of Class and Class I protostars (jShirlev et all 
120021 lYoung et alj|2003) . The isothermal dust tempera- 
ture is the temperature that, when used in an isothermal 
mass equation (e.g., Eq. [1]) yields the same mass as a 
detailed radiative transfer model including temperature 
gradients. There is a factor of 1.9 difference in mass be- 
tween assuming To = 10 K and 15 K. 

Figures [6] and [7] compare the peak and mean densities 
of starless and protostellar cores in each cloud. We use 
the peak column density Nh 2 , calculated from the peak 
1.1 mm flux density as a measure of the peak 

density: 



S, 



QbeamfJ>H2 m HK u B v (To) 



(2) 



with N(R 2 )/A V = 0.94 x 10 21 cm" 2 mag" 1 
(|Frerking et al.|[l982D . Here flbeam is the beam solid an- 
gle, niH is the mass of hydrogen, and ^h 2 — 2.8 is the 
mean molecular weight per H2 molecule. As discussed 
in Paper III, the 1.1 mm emission detected by Bolocam 
traces significantly higher column densities than other 
tracers such as the reddening of background stars. 

Mean particle densities, given in Tables H] and [5J are 
calculated within a fixed linear aperture of diameter 
10 4 AU: 

_ 3M le4 

nie4 " 4,r(ZW2) : W 1 j 

where D\ e ± is the aperture size (10 4 AU), M\ e ^ is the 
mass calculated from the 1.1 mm flux within that aper- 
ture, and fi — 2.33 is the mean molecular weight per 
particle. An aperture of 10 4 AU corresponds to approx- 
imately 40" in Perseus and Serpens, and 80" in Ophi- 
uchus. Note that n\ e ^ does not depend on the core size; 
while in some cases 10 AU may be considerably smaller 
or larger than the source FWHM, a fixed linear aperture 
is used here to mitigate the effects of linear resolution, 
which was found in Paper III to significantly bias the 
mean density calculated within the FWHM contour. 

The average peak column density (Nh 2 ) for starless 
cores in Perseus is 12 x 10 21 cm" 2 , while the average for 
protostellar cores is 50% higher, 19 x 10 21 cm" 2 , with a 
much wider dispersion (the significance of the Student's 
T-statistic is 0.03). Similarly, the typical mean density 
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Fig. 3. — Bolocam maps of Perseus, Serpens, and Ophiuchus, with the positions of starless and protostellar cores indicated. Identified 
cores have peak flux densities of at least 5<r, where a is the local rms noise level (on average a = 15 mjy beam -1 in Perseus, 10 mjy beam -1 
in Serpens, and 25 mjy beam -1 in Ophiuchus). Note that due to the large scales, individual structures are difficult to see. Regions of the 
maps with no detected sources have been trimmed for this figure. Starless and protostellar cores cluster together throughout each cloud, 
with both populations tending to congregate along filamentary cloud structures. There are a few exceptional regions, however, which are 
dominated by either starless or protostellar cores (e.g. the Bl Ridge, Serpens Cluster A). 

{ni e i) of starless cores (1.7x 10 5 cm~ 3 ) is a factor of three smaller than that of protostellar cores (5.1 x 10 5 cm~ 3 ; 
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Fig. 4. — Distributions of the angular deconvolved sizes of starless and protostellar cores in the three clouds. The size is given in units 
of the beam FWHM (0 m j,), and the mean of each distribution ± the dispersion in the sample is listed. Starless cores are larger on average 
than protostellar cores in Perseus, with a flattened distribution out to 3.5 m (,. In Serpens and Ophiuchus, however, there is little difference 
between the starless and protostellar distributions, and there are fewer very lar ge cores, particular ly in Serpens. The value of ddec/^mb can 
be used to infer the steepness of the source radial density profile (see text and I Young et ai]|2003ft . 
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Fig. 5. — Distributions of the axis ratios of starless and protostellar cores. Cores with axis ratios < 1.2 are considered round and those 
with axis ratios > 1.2 elongated, based on Monte Carlo simulations (Paper I). Only in Serpens, where starless cores tend to be slightly 
more elongated than protostellar cores, is there a distinguishable difference between the starless and protostellar populations. 

Ill, due in part to few sources near the detection limit 
and in part to the difference in calculating densities in a 
fixed aperture. 

Density distributions for starless and protostellar cores 
in Serpens are similar to those in Perseus. Peak Nh 2 
values are substantially smaller for starless ((Nh 2 ) = 
10 x 10 21 cm" 2 ) than for protostellar {(N H2 ) = 25 x 
10 21 cm -2 ; Student's T-test significance 0.04) cores, and 
form a much narrower distribution. Likewise, typical 
mean densities of starless cores (1.4 x 10 5 cm -3 ) are 
nearly three times smaller than those of protostellar cores 



Student's T-test significance 6 x 10~ 5 ). The large dif- 
ference in mean densities is due primarily to the sig- 
nificantly smaller sizes of protos tellar cores in Perseus. 
Recently, lJ0rgensen et al.l (|2007f ) found a similar result 
for Perseus, comparing SCUBA 850 /im cores with and 
without internal luminosity sources, as determined using 
Spitzer c2d data. Those authors concluded that cores 
with embedded YSOs (located within 15" of the core po- 
sition) have higher "concentrations" on average. Note 
that average ni e 4 values are much larger than the mini- 
mum detectable density of ~ 2 x 10 4 cm -3 from Paper 
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Fig. 7. — Distributions of th e m ean density ni e 4 of starless and protostellar cores, where the density is calculated in a fixed linear 
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lines, see Paper III). 
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(3.5 x 10 5 cuT 3 ; Student's T-test significance 0.03). In 
contrast to Perseus, however, mean density differences 
in Serpens are due entirely to the higher masses of pro- 
tostellar cores, as starless and protostellar cores have 
similar sizes in Serpens. As was the case for core sizes 
and shapes, there is essentially no difference between 
the peak or mean densities of starless and protostellar 
cores in Ophiuchus. Average peak Nh 2 values are sim- 
ilar (30 — 40 x 10 21 cm~ 2 ) and the two distributions 
have similar dispersions, and average ni e 4 values are 
2.2-2.8 x 10 5 cm" 3 . 

To summarize, peak column densities and mean densi- 
ties of starless cores in Perseus and Serpens are lower on 
average than for protostellar cores, whereas in Ophiuchus 
there is no significant difference between the starless and 
protostellar populations. The N^ 2 distributions for star- 
less cores are quite narrow in all three clouds, indicating 
a small range of column densities. This narrow distri- 
bution more likely represents an upper limit to Nh 2 for 
starless cores, as we are not sensitive to very large low 
mass cores (see § !3.3p . but in general cores with the high- 
est column densities tend to be protostellar. 

3.3. The Mass versus Size Distribution 

Figure [5] plots total core mass versus angular FWHM 
size for starless and protostellar cores in Perseus. The 
beam size and empirical 50% and 90% completeness lim- 
its as a function of size, derived from Monte Carlo simu- 
lations (Paper I), are indicated. Completeness limits are 
lower for protostellar cores by a factor of two due to the 
higher dust temperature assumed in the mass calculation 
(15 K) compared to starless cores (10 K). There is ap- 
proximately a factor of 1.2 increase in mass between the 
50% and 90% completeness curves, which holds true for 
the other clouds as well. 

In this diagram, starless cores seem to follow a con- 
stant surface density (M oc R 2 ) curve, consistent with 
the narrow distribution of peak column densities (Fig- 
ure [()]). Protostellar cores, in contrast, have a narrower 
range of sizes for a somewhat larger range of masses. 
This is a restatement of the results from § 13.11 and § 13.21 
protostellar cores in Perseus are smaller and have higher 
mean densities than starless cores. Note that because our 
completeness limits are similar to M oc i? 2 , the distribu- 
tion of starless cores in mass and size only implies that 
the upper envelope of cores follows a constant surface 
density curve. For example, there could be a population 
of large, low mass cores that we are unable to detect. 
Nevertheless, the two populations seem to fill different 
regions of the mass versus size parameter space. 

Examining Figure it is easy to imagine how proto- 
stellar cores in Perseus might have evolved from the cur- 
rent population of starless cores, by decreasing in size and 
increasing in density for a constant mass, until collapse 
and protostellar formation is triggered. Equivalently, the 
formation of a central protostar within a previously star- 
less core is associated with a decrease in core size and an 
increase in core density. 

Such a simple scenario is not consistent with the other 
clouds, however, as is evident in Figure [5] Although 
protostellar cores in Serpens have a small range in sizes 
for a large range of masses, as seen in Perseus, there is 
no population of large starless cores in Serpens. In fact, 
it is unclear how the relatively massive protostellar cores 



((M) PS = 2.1) in Serpens could have evolved from the 
current population of compact, low mass starless cores 
((M) SL = 0.9). It appears that Serpens has exhausted its 
reserve of starless cores with M > 2M Q ; unless new cores 
are formed from the lower density medium, future star 
formation in the cloud will result in stars of considerably 
lower mass than the current protostellar population. 

This deficit of starless cores at higher core masses may 
be the result of a mass dependence in the timescale for 
protostellar formation, with higher mass cor es forming 
protostars more quickly (|Hatchell et al.l 12008). In both 
Serpens and Perseus the most massive cores tend to be 
protostellar in nature, although the trend is much more 
extreme in Serpens. In Ophiuchus, where core masses 
are lower on average than in the other clouds, no such 
distinction between starless and protostellar cores is seen. 
There is essentially no difference between the starless and 
protostellar populations in Ophiuchus, suggesting very 
little core evolution after the formation of a protostar. 

3.4. Relationship to Cloud Column Density 

We use the cumulative fraction of starless and proto- 
stellar cores as a function of cloud Ay, shown in Fig- 
ure [lOl to quantify the relationship between dense cores 
and the surrounding cloud material. Visual extinction 
(Ay) is a measure of the cloud column density, and is de- 
rived based on the reddening of background 2MASS and 
IRAC sources, as described in Paper III and lHuard et al.l 
(2006). In Paper III, we found that 75% of 1.1 mm 
cores in Perseus, Serpens, and Ophiuchus are found at 
visual extinctions of Ay > 8 mag, Ay > 15 mag, and 
Ay > 20 — 23 mag, respectively. Although these values 
do not define a strict threshold, below these Ay levels the 
likelihood of finding an 1.1 mm core is very low. Here 
we investigate whether the relationship between dense 
cores and cloud column density is different for starless 
and protostellar cores. 

Figure [10] demonstrates that both starless and proto- 
stellar cores are found primarily at high cloud column 
densities (Ay > 6 mag): 75% of starless and proto- 
stellar cores in Perseus are located at Ay > 6.5 and 
9.5 mag, respectively, Ay > 8 and 10 mag in Serpens, 
and Ay > 19.5 and 25.5 mag in Ophiuchus. Note that 
these "threshold" values are significantly higher than the 
minimum cloud Ay which is approximately 2 mag in 
Perseus and Ophiuchus, and 6 mag in Serpens. Starless 
cores tend to be found at somewhat lower Ay than pro- 
tostellar cores; a two-sided KS test yields probabilities of 
approximately 1% that the starless and protostellar dis- 
tributions are drawn from the same parent distribution 
in each cloud. In Serpens and Ophiuchus there appear to 
be strict extinction thresholds for 1 . 1 mm cores at 7 and 
15 mag, respectively. A few cores in Perseus and Serpens 
lie outside the Ay map area, and may be associated with 
law-Ay material. 

As discussed in P aper III , an ex tinction threshold has 
been predicted by iMcKed (|1989t ) for photoionization- 
regulated star formation in magnetically supported 
clouds. In this model, core collapse and star formation 
will occur only in shielded regions of a molecular cloud 
where Ay > 4 - 8 mag. The fact that 75% of both pro- 
tostellar and starless cores are found above Ay ~ 6 mag 
in each cloud is consistent with this model; while it is 
certainly not the only explanation, magnetic fields may 
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Fig. 8. — Total mass versus angular FWHM size for starless and protostellar cores in Perseus. Dashed and dash-dot lines indicate 
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triggered. 
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there is no population of large starless cores in Serpens; it is not clear how the current population of relatively massive protostellar cores 
could have evolved from such compact, low mass starless cores. This discrepancy suggests that future star formation in the cloud may 
result in stars of considerably lower mass than the currently forming protostars. 
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Fig. 10. — Cumulative fraction of starless and protostellar cores as a function of cloud Ay, where the Ay is derived from the reddening 
of background stars in 2MASS and IRAC data. In all three clouds, the majority of cores are found at high cloud column densities 
(Ay > 7 mag). Dotted lines indicate the Ay at which > 75% of cores are found; 75% of starless and protostellar cores in Perseus are 
located at Ay > 6.5 and 9.5 mag, respectively. The equivalent values are Ay > 6 and 10 mag in Serpens, and Ay > 19.5 and 25.5 mag in 
Ophiuchus. There appears to be a strict extinction threshold in Serpens and Ophiuchus, with no cores found below Ay ~ 7 and 15 mag, 
respectively. 
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play a role in inhibiting collapse of cores, at least in the 
low column density regions of molecular clouds. 

3.5. Clustering 

Finally, we look at the spatial clustering of starless 
and protostellar cores. The spatial distributions of star- 
less and protostellar cores in each cloud are shown in 
Figure [3] Starless and protostellar cores appear to clus- 
ter together for the most part, often congregating along 
filamentary structures in the clouds. Both populations 
occur in the known clusters (e.g. NGC 1333, L 1688) and 
in a more distributed way across the clouds, with a few 
exceptions, such as in Perseus where there is a group of 
primarily starless cores (Bl ridge). 

As a more quantitative measure of clustering, we use 
the two-point correlation function (Figure [TT]): 



w ( r > = 7FTT ~ 1; 
H r (r) 



(4) 



where H s (r) is the number of core pairs with separa- 
tion between log(r) and log(r + dr), and H r (r) is similar 
but for a random distribution. The correlation function 
w(r) is a measure of the excess clustering as compared 
to a random distribution of sources. The upper panels 
of Figure [TT] plot w(r) as a function of source separation 
r, with the linear beam size and average source size indi- 
cated. The lower panels plot log(w), with a power law fit 
(w(r) oc r p ) for r larger than the average source size. In 
essence, the amplitude of w(r) is a measure of the mag- 
nitude of clustering, while the slope is a measure of how 
quickly clustering falls off on increasing scales. 

The amplitude of w(r) is higher for the protostellar 
samples in all three clouds, indicating that the degree 
of clustering is stronger on all spatial scales for proto- 
stellar cores. Visually, however, starless and protostellar 
cores tend to cluster in a similar way (Figure [3]), and 
this observation is supported by the similarity in the 
slope of w(r) for starless and protostellar cores. In Ser- 
pens, p = —1.4 ± 0.3 and —1.3 ± 0.6 for the starless 
and protostellar populations, respectively. The Ophi- 
uchus curves are quite noisy, but the best fit slopes 
(—1.3 ± 0.3 and —2.2 ± 0.4) are consistent within 2a. 
For Perseus the starless correlation function is shallower 
by 4a (p = -1.14 ± 0.13 and -1.5 ± 0.1). A shallower 
slope suggests that, while the amplitude of clustering is 
weaker for starless cores, it does not fall off as fast at 
larger spatial scales. 

A lower amplitude of clustering for starless cores as 
compared to protostellar cores in Perseus and Serpens is 
confirmed by the peak number of cores per square par- 
sec: 10 pc~ 2 and 16 pc -2 , respectively, for starless and 
protostellar cores in Perseus, and 4 pc -2 and 8 pc~ 2 in 
Serpens. These numbers suggest that clustering in the 
protostellar samples is a factor of 1.5-2 times stronger 
than in the starless samples. In Ophiuchus, however, the 
values of the starless and protostellar populations are 
identical (12 pc -2 ). 

There are at least three plausible reasons that cluster- 
ing might be stronger for protostellar cores, two envi- 
ronmental and one evolutionary. The difference may be 
an environmental effect if cores that are located in re- 
gions of the cloud with higher gas density are more likely 
to collapse to form protostars. In that case, more clus- 
tered sources would tend to be protostellar rather than 



starless. Similarly, if outflows or other protostellar activ- 
ity trigger the collapse of nearby cores, we would again 
expect more protostellar sources in clustered regions. If 
evolution plays a more important role, on the other hand, 
the spatial distribution of cores might evolve after pro- 
tostellar formation, for example as a result of dynamical 
effects. 

4. ARE THE STARLESS CORES REALLY PRESTELLAR? 

For any discussion of the mass distribution or lifetime 
of prestellar cores, it is important to determine whether 
or not our starless cores are likely to be truly prestellar 
(i.e. will form one or more stars at some point in the 
future). To this end, we investigate the dynamical state 
of starless cores by estimating the virial masses of cores 
in Perseus. 

While it is possible that some fraction of the observed 
starless cores are transient or stable structures that will 
never form stars, the starless cores detected by our Bolo- 
cam surveys have high mean particle densities (typical 
n le4 ^1 — 3 x 10 5 cm" 3 ; £13.21), mak i ng them likely 
to be prestellar (|Di Francesco et al.ll2007l ; lKeto fc Casellil 
2008). Comparison of our data to molecular line ob- 
servations is a more robust method of determining if 
cores are gravitationally bound; line- widths provide a di- 
rect estimate of the internal energy of cores, while the 
1.1 mm dust masses provide an estimate of the poten- 
tial energy. An ammonia (NH3) (J,K)=(1,1) and (2,2) 
survey of dense cores in Perseus that includes all of the 
Bolocam- identified 1.1 mm cores h as recently been com- 
pleted bv lRosolowskv et al.1 (|2008f ) at the GBT as part of 
the COMPLETE project. 10 These observations provide 
a much-needed measure of the temperature and internal 
motions of dense millimeter cores. 

Figure [12] plots the ratio of the dust mass (Md ust , 
from Tables 0] and [5]) to the virial mass (M v ; r ) as 
a function of Mdust f° r the 70 cores in Perseus that 
have well-determined lin e widths and temperatures from 
(|Rosolowskv et al.1 ^008) NH 3 survey. Although all Bolo- 
cam cores were detected in ammonia, not every target 
had a significant NH3 (2,2) detection or sufficient signal- 
to-noise to measure the line width in the presence of the 
hyperfine structure. We present the data relevant to Fig- 
ure [TJ in Table H 

We estimate the virial mass of the cores as 



aG ' 



(5) 



where a accounts for the central concentration of the 
core, R is the radius of the core, and <7h2 is the line width 
of the molecular gas inferred from the ammonia line 
width. Note that Mdust/-M vir = 1/a for a = 1, where a 
is the virial parameter introduced bv lBertoldi fc McKed 
(1992). For simplicity, we adopt a power-law density 
profile, for which a = (1 — p/3)(l — 2p/5)~ 1 . We report 
values in Table [2] for a,p = 1.5 profile, giving a = 5/4; for 
p = 2 (SIS) the dynamical masses would be 25% smaller. 
The total line width of H2 includes thermal and turbulent 
components, and is calculated as 

2 2 kT K kT K 



'H2 



NH3 



\lm\ 



2mn 



(6) 



10 "The COMPLETE Survey of Star For ming Regions"; 
http://cfarwww.harvard.edu/COMPLETE/] HGoodman et al.1 
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Fig. 11. — Two-point spatial correlation function w(r), plotted as a function of source separation r, for starless and protostellar cores 
in Perseus, Serpens, and Ophiuchus (upper panels). Lower panels plot log(«>), with the best fitting power law slope for each distribution. 
Resolution limits and average source sizes are indicated. Power law slopes are fit for r larger than the average source size in each cloud. The 
amplitude of the protostellar w(r) is consistently higher than the starless w(r), indicating that clustering of protostellar cores is stronger 
on all spatial scales. Slopes for the starless and protostellar populations are similar in each cloud, however, so although the magnitude of 
clustering differs, the fundamental nature of the clustering does not. 



The radius is estimated from the 1.1 mm data as 
R = 9 dec /2, or the HWHM. Although the effective ra- 
dius R appropriate for the virial theorem is not exactly 
equivalent to the gaussian-fit HWHM, we do not expect 
the correction to be large. If the FWHM is a good repre- 
sentation of v<Bhi2 times the radial dispersion one would 
measure in the x-y plane for a power law density profile, 
then we expect 

R = (WV(5 - p)/(3 - p) % /37Vl/(81n2) (7) 
or R ~ 0.56 Odec for p — 1.5, only a 10% correction. 

Dotted and solid lines in Figure [T2] indicate regions in 
the plot for which cores should be virialized and self- 
gravitating, respectively. Note here that Md us t = M V [ T 
implies 2K = — U, where K is the kinetic and U the 
gravitational energy, while the self-gravitating limit is 
defined by K = —U, or Md us t/M v i r = 0.5. There are 
two important points to note here. First, nearly all of 
the starless cores lie close to or above the self-gravitating 
line, indicating that they are likely to be gravitationally 
bound. This statement must be modified by the impor- 
tant caveat that the power law profile assumed and the 
method for measuring the dust mass introduce significant 
systematics that could shift Mdust/ M v - lT by more than a 
factor of two. 

Perhaps more convincingly, there is very little differ- 
ence between starless and protostellar cores in this plot; 
even though several of the starless cores lie below the self- 



gravitating line, this regions is populated by protostellar 
cores as well, which are by definition capable of forming 
stars. F urthermore, recent estimates of the dust opacity 
«i.i mm (|Shirlev et al.ll2007l see ? I3.2[) may increase Md us t 
by 1.3, making the starless cores more bound. Based on 
this discussion, we assume from here on that all starless 
cores in our 1.1 mm samples are true prestellar cores. 

5. THE PRESTELLAR CORE MASS DISTRIBUTION AND 
THE IMF 

One very important measure of the initial conditions 
of star formation is the prestellar core mass distribution 
(CMD). In particular, comparing the prestellar CMD to 
the stellar initial mass function (IMF) provides insight 
into how the final masses of stars are determined. There 
are a number of processes that may (jointly) dictate what 
the final mass of forming star will be. We focus here on 
two simple cases, assuming that only one drives the shape 
of the IMF. 

If stellar masses are determined by the initial fragmen- 
tation into cores, i.e., the final star or binary mass is 
always a fixed percentage of the original core mass, then 
the shape of the emergent stellar IMF should closel y trace 
that of the prestellar CMD (e.g-. lMvers et aljl 998). This 
might be expected in crowded regions where the mass 
reservoir is limited to a protostar's nascent core. If, on 
the other hand, stellar masses are d etermined by compet- 
itive accretion (Bo nnell et al.|[200l[ ). or by the protostars 
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Dynamical Properties of Perseus Cores 
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Note. — Kinetic temperatures {Tk ) and lincwidths (<th2 ) are 
derived from the NH3 observations of Perseus cores completed by 
IRosolowsky et all fl 20081) . The virial mass Af v ; r is calculated accord- 
ing to Eq. [5] with the effective radius R estimated as the deconvolved 
HWHM of the 1.1 mm identified core (#d ec /2). Af^ust is the mass 
calculated from the total 1.1 mm flux (Tables [41 and [5j . 
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Fig. 12. — Ratio of the total mass derive from dust (M^ust) 
to the virial mass (M v i r ), calculated using the temperature and 
velocity dispersion derived from GBT NH3 (1,1) and (2,2) obser- 
vations (Rosolowskv et al. 2008), for starless and protostellar cores 
in Perseus. Only those with well-determined NH3 line- widths are 
shown. Solid and dashed lines indicate the minimum Mdust/M V i T 
for which cores should be in virial equilibrium (dashed) or self- 
gravitating (solid). 

themselves through feedback mechanisms (e.g., outflows 
and winds; IShu et al.|[l987[) . we would not expect the 
emer gent IMF to reflect the original core mass distribu- 
tion (I Adams k. Fatuzzo|[l99l . 

5.1. The Prestellar CMD 

We combine the prestellar core samples from all three 
clouds, 108 cores in total, and assume To = 10 K to 
calculate masses. As noted above in § 01 our assump- 
tion that the majority of starless cores in our sample are 
truly prestellar is supported by a comparison of the dust 
mass to the dynamical mass from NH3 observations. The 
resulting prestellar CMD is shown in Figure [T51 

The 50% completeness limit shown (M ~ 0.8 M Q ; dot- 
ted line) is estimated based on the empirical 50% com- 
pleteness curve in Figure [5J Our completeness depends 
on the size of a given source; as the average core size 
in Perseus is 68", we take the 50% completeness limit 
for a 70" FWHM source (M - 0.8 M ) to be the aver- 
age 50% completeness limit of the sample. The equiv- 
alent 90% completeness limit is a factor of 1.2 higher, 
M ~ 0.93 Mq. 50% completeness limits for average-sized 
sources are lower in Serpens (0.6 Mq) and Ophiuchus 
(0.5 Mq), but as more than half of the total population 
of prestellar cores are in Perseus, we take 0.8 M for the 
entire prestellar sample. 

We fit a power law (dN/dM oc M a ) to the CMD 
for M > 0.8 Mq, finding a slope of a = —2.3, with 
a reduced chi-squared of x 2 = 1.9. The best fit slope 
depends somewhat on the histogram binning, ranging 
from 2.0 to 2.6 for bin widths of 0.1 to 0.3 Mq, so we 
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Fig. 13. — Combined prestellar core mass distribution (CMD), with power law and lognormal fits. The prestellar sample is composed of 
all starless cores from Perseus, Serpens, and Ophiuchus, and the 50% mass completeness limit (dotted line) is defined by the completeness 
limit for average-sized cores in Perseus. Recent measurements of the ste l lar IM F for M > 0.5M p) (a = —2.3 to —2.8) are similar to the 
best-fit CMD power law slope (a = —2.3 ±0.4). IMF fits from Chabricr (2005) (lognormal) and Kroupa (2002) (three-component power 
law) are shown as thick gray lines for reference. The shad ed h istogram indicates how the mass distribution changes if a small fraction of 
cores (15%, based on the six "unbound" cores from Figure fT2l with M < IMq are excluded from the prestellar sample. 



assign a total uncertainty of 0.4 to our best fit slope, 
when also taking into account formal fitting errors. We 
fit a lognormal distribution to M > 0.3 M©, finding a 
best-fit width a = 0.30 ± 0.03 and characteristic mass 
Mo = 1.0 ± 0.1 Mg. Although the lognormal function is 
quite a good fit (x = 0.5), the reliability of the turnover 
in the prestellar CMD is highly questionable given that 
the completeness limit in Perseus coincides closely with 
the turnover mass. The prestellar CMD can also be fit by 
a broken power law with a = — 4.3± 1.1 for M > 2.5 Mq 
and a = -1.7 ±0.3 for M < 2.5 M , although the un- 
certainties are large. 

Given that source detection is based on peak intensity, 
we may be incomplete to sources with very large sizes 
and low surface density even in the higher mass bins 
(M > O.8M ). Completeness varies with size similarly 
to M oc i? 2 , thus the fraction of (possibly) missed sources 
should decrease with increasing mass. 11 The effect of 
missing such low surface brightness cores, if they exist 
and could be considered prestellar, would be to flatten 

11 Unless M oc R 2 intrinsically for starless cores, in which case 
a constant incompleteness fraction would apply over all mass bins. 



the CMD slightly (i.e. the true slope would be steeper 
than the observed slope). Instrumental selection effects 
are discussed further in Paper I. 

To be completely consistent, we should exclude 
the "unbound" cores from Figure [TJ] (those with 
M dust /M vir < 0.5) from our prestellar CMD. This rep- 
resents 6 out of the 40 cores that have measured virial 
masses in Perseus, all 6 of which have Md us t < 1M Q . 
We do not have virial masses for cores in Serpens or 
Ophiuchus, but we can randomly remove a similar frac- 
tion of sources with M < 1M Q from each cloud sam- 
ple (2 sources from Serpens, 4 from Ophiuchus, and an 
additional 4 from Perseus). The shaded histogram in 
Figure fl~3l indicates how the mass distribution is altered 
when these 16 "unbound" cores are excluded from the 
sample. Our derived CMD slope is not affected, as nearly 
all of the starless cores below the "gravitationally bound" 
line in Figure [12] have masses below our completeness 
limit, and even at low masses the CMD is not signifi- 
cantly changed. 

There may also be some concern over the use of a sin- 
gle dust temperature Tr> = 10 K for all cores. To test 
the validity of this assumption, we use the kinetic tern- 
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Fig. 14. — Effect on the presteffar CMD in Perseus of using NH3- 
derived core kinetic temperatures (shaded histogram) rather than 
assuming a constant Tjj = 10K for all cores (black histogram). 
Only starless cores in Perseus are included here. The best-fitting 
power-law slope for M > O.8M0 is similar for the two histograms: 
a = -2.35 ± 0.33 and a = -2.33 ± 0.28 for the T D = 10 K 
and NH3 temperature curves, respectively. Although we do not 
have kinetic temperatures for cores in Serpens or Ophiuchus, we 
do not expect temperature effects to dramatically alter the derived 
prestellar CMD slope. 

perature s (Tk) derived from th e GBT NH3 survey of 
Perseus (jRosolowskv et all 120081 Schnee et al. 2008, in 
prep) to compute core masses, assuming that the dust 
and gas are well coupled (i.e., Tr> = Tk)- Figure ITU 
shows the CMD of prestellar cores in Perseus, both for 
a single To = 10 K, and for masses calculated using the 
NH3 kinetic temperatures for each core. There is some 
change to the shape of the CMD at intermediate masses, 
but the best-fitting slope for M > O.8M (a = -2.3) is 
unchanged. In fact, the deviation from a power law is 
smaller when using the kinetic temperatures than when 
using Td = 10 K (x 2 — 1.3 and 2.5, respectively). 

The median Tk of prestellar cores in Perseus is 10.8 K, 
quite close to our adopted Tp = 10 K, and the small 
overall shift in masses (a factor of 1.1) corresponding to 
an 0.8 K temperature difference would not affect the de- 
rived CMD slope. The dispersion in kinetic temperatures 
in Perseus is ±2.4 K, or ±0.4M Q for a 1M© core, and the 
tail of the distribution extends to Tk > 15 K (Schnee et 
al. 2008, in prep). We do not have temperature infor- 
mation for cores in Serpens or Ophiuchus; if the median 
temperature were to vary from cloud to cloud by more 
than a few K, the shape of the combined CMD could be 
significantly altered. Perhaps a more serious issue is that 
prestellar cores are not in reality isothermal, but taking 
into account variations of the temperature with radius 
requires radiative transfer modeling of each source. 

5.2. Comparison to the IMF 
The shape of the local IMF is still uncertain (jScalol 



2005) , but recent work has found evidence for a slope of 
a = —2.3 to —2.8 for stellar masses M > 1 M©, simi 



lar to the slope we measure for the c ombined prestellar 
CMD (a = -2.3 ± 0.4). For example. [Reid et al.l (|2002f ) 
find a = —2.5 above 0.6 M© , and a — —2.8 above 
1 M©. ISchroder fe Pagell (|2003f ) suggest a = -2.7 for 
l.K M < 1.6 M© and a = -3.1 for 1.6 < M < 4 M Q . 
For reference, t he Salpete r IMF has a slope of a = —2.35 
(|Salpeterl lT955). and the IScalol (1986) slope for sources 
with mass M > 1 M© is a ~ —2.7. At lower masses, 
the IMF flatt ens, and may b e characterized by a lognor- 
mal function. iKroupal ()2002f ) suggests a three-component 
power law for the average single-star IMF: a = —2.3 for 
0.5 < M < 1 M©, a = -1.3 for 0.08 < M < 5 M^ 
and a = -0.3 for 0.01 < M < 0.08 M©. iChabrierl (p005h 
finds that a lognormal distribution with a = 0.55 and 
M = 0.25 Mp) is a go od fit for M < 1 M©. 

T he IKroupal (I2002D three-component power law and 
the lChabrieiT (i2005) lognormal IMFs are shown as thick 
gray lines in Figure [T3l The width of the lChabrieil (|2005l ) 
lognormal (a = 0.55) is somewhat larger than the width 
of the prestellar CMD best-fit lognormal (a = 0.3), as 
expected if we are incomplete at lower masses, and the 
IMF characteristic mass (Mo — 0.25 M©) is a factor of 
four lower than that of the CMD (M = 1.0 M©). A 
lower characteristic mass for the IMF is expected if some 
fraction of the core mass is lost in the star-formation 
process. Power law fits for M > 1 M© appear to be 
quite similar for the CMD [a = -2.3 ± 0.4) and IMF 
(a = —2.3 to —2.8), however. 

Although we cannot rule out the importance of feed- 
back and other local processes in determining the shape 
of the IMF, the fact that the prestellar CMD and the 
local IMF have similar shapes supports a growing body 
of evidence that the final masses of stars are determined 
during the core formation process. If this is the case, it 
is tempting to relate both the prestellar CMD and the 
IMF to the CMD created from turbulent simulations. 
Numerical models of turbulent clouds have had some suc- 
cess in reproducing the general shape o f the IMF (e.g., 
iPadoan fc Nordlundl[2002l : iLi et al.ll2004D . but the link is 
not firmly established. For example, we found in Paper 
III that the predicted de pendence of CMD shape on the 
turbulent mach number (|Ballesteros-Paredes et al.ll2006l : 
IPadoan fc Nord lund 2002]) does not agree with observa- 
tions. 

Evidence for a direct link between the CMD and the 
IMF has been found previousl y based on dust emis- 
sion surveys of small regions (|Testi fc Sargentl 119981 : 
iMotte et al.lll998[). as well as mol ecular line observations 
of den se cores (jOnishi et aLll2002f h Recently. I Alves et al.1 
(|2007h found evidence for flattening at low masses in 
the CMD of the Pipe Nebula, as traced by dust extinc- 
tion toward background stars. Those authors interpret 
the similarity between the Pipe Nebula CM D and the 
Trapezium cluster IMF (jMuench et al.ll2002] ) as evidence 
that the stellar IMF is a direct product of the CMD, 
with a uniform core-to-star efficien cy of 30% ± 10% . 
Although the measured masses of lAlves et al.1 J2007|) 
are somewhat less uncertain than ours because they do 
not need to assume a dust opacity or temperature, the 
mean particle densities of cores traced by dust extinc- 
tion (n~ 5x 10 3 — 2 x 10 4 cm~ 3 ) are considerable lower 
than the mean densities of cores traced by our Bolocam 
1.1 mm surveys (n ~ 2 x 10 4 — 10 6 cm" 3 ), and they 
may never form stars. In fact, recent C 18 and NH3 
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observations suggest that the dust extinction sources in 
the Pipe Nebula are p ressure confined, gravitationally 
unbo und starless cores (|Lada et alj l2008t iMuench et al.l 
12007ft . In Orion. INutter fc Ward- Thompson! (|2007t ) find 
a turnover in the CMD of starless SCUBA 850/im cores 
at ~ 1.3M©. T hose authors re late this turnover to a 
down-turn in the iKroupal (|2002f ) IMF at - 0.1M Q , and 
infer a much lower core-to-star efficiency of ~ 6%. 

If the prestellar CMD does have a one-to-one relation- 
ship with the stellar IMF, then the ratio of turnover 
masses is a measure of the core collapse efficiency, or 
the fraction of original core mass that ends up in the fi- 
nal star: f cB = M^ F /Af^ D . Here M TO is the mass 
where the dN/dM distribution, which rises with decreas- 
ing mass, flattens out and begins to fall with further de- 
creasing mass. Equivalently, M^ D /M^ F = 1 — / off 
is the fraction of core mass lost in the star formation 
process. 

Our limited completeness dictates that we can only 
measure a lower limit to f c g . If there is a true turnover 
in the prestellar CMD, it must occur below our complete- 
ness limit, at Mg^ D < 1.0 M Q . The system IMF, i.e., 
treating binaries and multiple systems as single rather 
than multiple objects peaks at M^ F ~ 0.2 - 0.3 M Q 
(e.g.. IChabrien 120051 iLuhman et alll2003ft . This system 
IMF is appropriate for comparison to our CMDs, as we 
would not resolve such multiple systems even if they 
form from distinct cores. For an IMF turnover mass of 
0.25 M Q , M^ d < 1.0 M implies that at least 25% 
of the initial core mass is accreted onto the final star 
or stellar system. Characteristic masses associated with 
lognormal fits to both the IMF and CMD imply a similar 
ratio of f cS > A^ MF /A/ CMD ~ 0.25/1.0 ~ 0.25. 

Our conclusion that the core-to-star efficiency is 
at least 25% is c onsistent with the value found by 
lAlves et al.l (|2007l ) for the Pipe Nebula (30%), and 
with predicted efficiencies of 25% — 75% from re- 
cent analytic models of bipolar protostellar outflows 
dMatzner fc McKedl2005h . 

5.3. The Effect of "Multiple" Sources 

As noted in § [3J 17% — 55% of our protostellar cores 
are associated with more than one (two to three) em- 
bedded protostars. We can expect that a similar frac- 
tion of prestellar cores will form a resolved binary or 
multiple system, making a direct mapping between the 
CMD and IMF difficult to justify unless the resolutions 
are matched. Given the relatively low resolution of the 
Bolocam data (30"), some cores must result in wide- 
separation multiple stellar systems that would not be 
considered single objects in the system IMF. Further- 
more, the system IMF likely evolves over time due to 
dynamical effects such as the decay of hi gh-order multi- 
ples, ejections, and close interactions (see lGoodwin et all 
20081 for more discussion of the CMD to IMF mapping) . 

The effect of these multiple sources will likely be to 
steepen the CMD, as higher mass cores are divided into 
multiple lower mass sources. This is consistent with the 
fact that the prestellar CMD slope is at the low end of 
the IMF slope range. It is important to keep in mind 
that we would expect only ~ 25 — 30% of cores in the 
combined prestellar CMD to be "multiple" sources, so 
they should not dominate the CMD sl ope. While we 
know that most stars occur in binaries (|Duchene et al.l 



120071 and references therein), and thus that there must 
be some unresolved binaries in the Spitzer data, these 
compact systems will most likely be unresolved in the 
system IMF as well, so they do not affect our comparison. 

5.4. The Protostellar Core Mass Distribution 

For comparison, we show in Figure [T5] the mass dis- 
tribution of the combined protostellar core sample from 
all three clouds, assuming Try — 15 K. The protostellar 
CMD is considerably wider and flatter than the prestel- 
lar CMD, and extends to higher masses. The best-fitting 
power law slope (a — —1.8; x 2 = 1.7) is shallower than 
for the prestellar sample, and the best-fitting lognormal 
distribution (a = 0.51 ± 0.07, M = 0.8 ± 0.2; x 2 = 1.3) 
is wider by nearly a factor of two. The protostellar CMD 
must be interpreted in the context of a constant assumed 
To = 15 K; an intrinsic spread of temperatures would 
tend to flatten the mass distribution even more, but 15 K 
has been determined to be the best single Tp based on 
many detailed radiative transfer models (|Shirlev et al.l 
l2lMlYmmg et al.ll2003ft . 

A protostellar CMD that extends to lower masses than 
the prestellar CMD is not unexpected, given that some 
fraction of the core mass has already been accreted on 
to the central source for protostellar cores. An extension 
to higher masses in the protostellar CMD as compared 
to the prestellar CMD, however, can only be explained 
if (a) we are underestimating the dust temperature of 
protostellar sources or underestimating the total mass of 
prestellar cores, (b) the current populations of prestellar 
cores in these clouds have lower mass than the generation 
of cores that have already formed protostars, or (c) the 
highest mass prestellar cores collapse to form protostars 
v ery quickly. 

IHatchell et all £2008) examine the relationship be- 
tween starless and protostellar mass distributions in 
Perseus based on SCUBA data, finding that a sim- 
ple mass-dependent evolutionary model in which higher 
mass cores form protostars on a shorter timescale than 
lower mass cores can qualitatively explain the observed 
steeper prestellar CMD slope, and extension of the pro- 
tostellar CMD to higher masses. Another step toward 
resolving the reasons for these differences may be careful 
radiative transfer models of each prestellar and proto- 
stellar source to replace the assumption of a constant 
T D . 

6. LIFETIME OF THE PRESTELLAR PHASE 

We can use the relative number of prestellar cores 
and embedded protostars to estimate the lifetime of the 
prestellar core phase, which is essential for understand- 
ing the physical processes that govern the formation, sup- 
port, and collapse of star-forming cores. Several assump- 
tions are required for this calculation. We assume that 
all of the starless cores in our sample are in fact prestel- 
lar, and will eventually collapse to form stars or brown 
dwarfs (see @. As we ultimately calibrat e our lifetimes 
based on the Cla ss II phase (2 x 10^ yr; iKenvon et alJ 
119901: ICieza et al.ll2007t ISpezzi et al.ll2008h . star forma- 
tion must have been steady in time for at least the last 
2 Myr. In addition, we must assume that there is no 
significant evolutionary dependence on source mass. Al- 
though this kind of analysis relies on a number of as- 
sumptions, it has the advantage of being quite simple, 
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Fig. 15. — Combined protostellar mass distribution, with power law and lognormal fits. The protostellar CMD includes all protostellar 
cores from Perseus, Serpens, and Ophiuchus, and the 50% mass completeness limit is defined by the completeness to average sized 
protostellar c ores in Perseus. The best-fit power law (a = —1.8) is shallower, and the distribution as a whole is wider, than the prestellar 
CMD (Figure |13> . This difference is expected if the protostellar CMD evolved from the prestellar CMD. The fact that the protostellar 
distribution extends to higher masses than the prestellar CMD is more difficult to explain, however, and may suggest that higher mass 
prestellar cores are relatively short lived. 



and it does not require on an a priori accretion rate or 
star-formation model. 

Table [3] lists the ratio of the number of starless cores 
(AT SL ) to embedded protostars (N cmh = ^ciassO+iVciassi) 
for all three clouds, and the starless core lifetime (<sl) 
derived from that ratio: t$L — t cm b(NsL/N cm b). The 
number of embedded protostars in each cloud is derived 
in a companion paper (Enoch et al. 2008, in prep), in 
which the SED derived from c2d Spitzer photometry and 
1.1 mm fluxes is used to classify protostars based on their 
bolometric temperature, Tboi- Class (Ttoi < 70 K) and 
Class I (70 < Tboi < 650 K). In addition to the T M 
classification, all embedded protostars are required to be 
detected at 1.1 mm, so that we are complete to embedded 
sources with M env > 0.1M Q (Enoch et al. 2008, in prep). 

The ratio of the number of starless cores to the number 
of embedded protostars is NsL/N em b = 1.0 in Perseus, 
0.4 in Serpens, and 0.8 in Ophiuchus. For an embedded 
protostar phase that lasts i e mb = 5.4 x 10 5 yr (Evans et 
al. 2008, in prep), the observed ratios imply prestellar 
core lifetimes of 5 x 10 5 yr in Perseus, 2 x 10 5 yr in 
Serpens, and 5 x 10 5 yr in Ophiuchus. Note that we 
also found approximately equal numbers of protostellar 
and starless cores in all three clouds (Table [T|) , further 
confirmation that the lifetime of starless cores is similar 
to that of the embedded protostellar phase. 12 Therefore, 

12 The N SL /N cmb ratios in Table |3] differ from N SL /N PS in Ta- 
ble [T] due to the presence of multiple embedded protostars in some 



the dense starless cores we are sensitive to last for 2 — 5 x 
10 5 yr in all three clouds. The differences in NsL/N cm b 
from cloud to cloud may be environmental. The ratio is 
lowest in Serpens, which has the highest mass cores on 
average. As noted in $5Al higher mass cores tend to be 
protostellar, possibly because they collapse on a shorter 
timescale than lower mass cores. 

Taking all three clouds together yields a prestellar core 
lifetime of 4.5 x 10 5 yr. The uncertainty in the measured 
number of prestellar cores is approximately ±10, based 
on the range in the number of starless cores for differ- 
ent identification criteria (see § 12. 3[) . Given a similar 
uncertainty in the number of embedded protostars (see 
Enoch et al. 2008, in prep.), this corresponds to an uncer- 
tainty in the lifetime of prestellar cores of approximately 
0.8 x 10 5 yr. 

Published measurements of the prestellar core life- 
time vary by two orders of magnitude, from a few times 
10 5 yr to 10 7 yr (I Ward-Thompson et all 120071 ). For ex- 
ample, iLee fc Mversl (|1999() calculate a core lifetime of 
6 x 10 5 yr for optically selec ted cores with mean densi- 
ties of 6 — 8 x 10 3 cm -3 , while [Jessop fc Ward-T hompson 
(2000) find a lifetime of 10 7 yr for low density cores de- 
tected from column density maps based on IRAS far- 

protostellar cores ([[3} and to the fact that some embedded pro- 
tostars are "band-filled" at 1.1 mm. Band-filled sources are not 
associated with a distinct core, but appear to be associated with 
millimeter emission that is either extended or below the 5<r detec- 
tion limit. 
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TABLE 3 

Lifetime of the prestellar core phase 



Cloud 


N SL /N cmb 


*SL 

(yr) 


<"le4>g L 
(cm -3 ) 


in 
(yr) 


*SL/*ff 


Perseus 


67/66 = 1.0 


5 x 10 5 


1.4 X 10 5 


1.4 X 10 5 


3.9 


Serpens 


15/35 = 0.4 


2 x 10 5 


1.2 X 10 5 


1.5 x 10 5 


1.5 


Ophiuchus 


26/28 = 0.9 


5 x 10 5 


1.6 x 10 5 


1.3 x 10 5 


3.9 


Combined sample 


108/129 = 0.8 


4.5 x 10 5 


1.7 x 10 5 


1.4 x 10 5 


3.2 


"low density" 


78/129 = 0.6 


3.3 x 10 5 


1.3 x 10 5 


1.5 x 10 5 


2.2 


"high density" 


30/129 = 0.2 


1.3 x 10 5 


2.8 x 10 5 


1.0 x 10 5 


1.3 



Note. — N cm b refers to the total number of embedded protostars, i.e., those in the Class and Class I phases: N cm b — Nci asa Q + Nciassl- The 
number of embedded protostars is derived in a companion paper (Enoch ct al. 2008, in prep). The prestellar core lifetime f-sL is derived from the 
number of starless cores: tsL — t C mb(^SL/^ cm b), for an embedded phase lifetime t om b — 5.4 x 10 5 yr (Evans et al. 2008, in prep). Core mean 
densities ni e 4 arc calculated in a fixed linear aperture of diameter 10 4 AU, and the free-fall timescalc iff of starless cores is derived from the typical 
mean density (ni e 4) SL of the starless samples using Eq. {8j. The "Combined sample" includes sources from all three clouds; this combined sample 
is divided into "low density" and "high density" bins representing cores with n± e 4 < 2 X 10 5 cm -3 and n± e 4 > 2 X 10 5 cm -3 , respectively. 



infrared observations. Ou r results are similar t o re- 
cent findings in Pers eus by iJorgensen et al.l (|2007( ) and 
lHatchell et al.l (|2007l ). both of whom find approximately 
equal lifetimes for the starless and embedded protostellar 
phases by comparing SCUBA 850 /im maps with Spitzer 
c2d data. The approximate equality between Nsl and 
./Vps has been fou nd by a number of stu dies, and was 
noted e arly on by iBeichman et alj (1986) for the NH3 
cores of iMvers fe Benson! (|1983) . Our results arc also 
consiste nt within a factor of two with the lifetime de- 
rived bv IVisser et a l. (2002) for a sample of Lynds dark 
clouds observed with SCUBA. 

The average mean density (calculated in a fixed linear 
aperture of 10 4 AU, see § I3.2[) of the starless core samples, 
("■ie4)gL! an d corresponding free-fall timescale, tg, are 
also given in Table [3] The free-fall time is the timescale 
on which starless cores will collapse in the absence of 
internal support, and is calculated from the mean particle 
density (|Spitzerlll97cl : 



iff 



3tt 
32Gp 



3tt 



32G(n) pLpiUH ' 



(8) 



where run is the mass of hydrogen and \x p = 2.33 is the 
mean molecular weight per particle. Mean densities are 
similar in all three clouds, (ni e 4) SL — 1 — 2 x 10 5 cm -3 , 
with corresponding free-fall times oftff = 1.3 — 1.5 x 
10 5 yr. The final column in Table [3] gives the ratio of 
the measured starless core lifetime to the average free- 
fall timescale in each cloud: tsh/tg = 1.5 — 3.9. Thus 
prestellar cores last for only a few free-fall times in all 
three clouds. 

Such a short prestellar core lifetime argues for a dy- 
namic, rather than quasi-static, core evolutionary sce- 
nario. One classical quasi-static model is that of magnet- 
ically dominated star formation, in which the evolution of 
highl y sub-critical cor es is moderated by ambipolar diffu- 
sion (|Shu et aT1ll987f ). In this paradigm, prestellar cores 
should have lifetimes similar to the ambipolar diffusion 
timescale, or £ad ~7x 10 6 yr for typical ionization lev- 
els in l ow-m ass star forming regions (e.g., lNakanolll998l : 
lEvansl[l999h . more than an order of magnitude longer 
than our results. 

It is important to emphasize here that our Bolocam 
surveys are sensitive to to cores with relatively high 



mean density (n > 2 — 3 x 10 4 cm -3 ; Paper III). 13 
Thus we may be sampling only the densest end stage 
in a longer core evolutionary picture, in which case a 
magnetic field dominated scenario could still be appl i- 
cable at early t imes (e.g., iTassis fc Mouschoviasll2004l) . 
ICiolek fc Basul (|2001[) also note that t^D can be as short 
as a few free-fall times for marginally sub-critical cores. 

6.1. Density Dependence 

Interestingly, if we divide our sample into two density 
bins, ni e 4 < 2 x 10 5 cm -3 and n ie 4 > 2 x 10 5 cm" 3 , we 
find a longer lifetime for lower density cores: 3.3 x 10 5 
yr versus 1.3 x 10 5 yr for the higher density sample (see 
Table [3|) . This result suggests that the prestellar core 
lifetime becomes shorter as cores become more centrally 
condensed, at a faster rate than n -0 5 (i.e., that expected 
based on the dependence of tg on n, see Eqj8]). The ob- 
served trend, isL oc n" 1 - 2 , is closer to the n~°- 8 5 dep en- 
dence suggested by Uessop fc Ward- Thompson! |2000) . 



6.2. Additional Uncertainties 

The prestellar lifetime depends on the ratio with proto- 
stellar cores, but our sensitivities to prestellar and proto- 
stellar sources are not equal due to the different assumed 
dust temperatures: the protostellar core mass complete- 
ness limit is lower by approximately a factor of two. Of 
the 55 protostellar cores in Perseus, 13 fall between the 
protostellar and starless completeness limits (similarly 
3/20 in Serpens, and 5/17 in Ophiuchus), suggesting that 
the prestellar core lifetime should perhaps be increased 
by a factor of 1.3. 

We must also consider the possible bias introduced by 
using the number of embedded protostars, as detected 
with the higher resolution Spitzer data, rather than the 
number of Bolocam-detected protostellar cores, to calcu- 
late the prestellar lifetime. N cn ±,/Npg = 1.2, 1.8, and 
1.6 in Per, Ser, and Oph, respectively, as can be seen by 
comparing Tables [T] and [3l If the starless cores have a 
similar " multiplicity" fraction that we missing due to the 
30" resolution of Bolocam, then the prestellar lifetimes 

13 Note that this limiting density takes into acc ount the 50% 
completeness limit to starless cores as discussed in £13.31 and i|5,ll 
it is the typical density found when one calculates a mean density 
along the 50% completeness curve as a function of size in Figures [8] 
and [9] 
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would be increased by these factors. If cores fragment 
at later times, or if multiple protostars collapse from a 
single core, our lifetimes would not be affected. 

On the other hand, if some fraction of the starless core 
sample are not prestellar but unbound starless cores, our 
calculated lifetime of 4.5 x 10 5 yr would be an overesti- 
mate of the true prestellar core lifetime. Extrapolating 
from the six cores with Md U st/-Mvir < 0.5 in Figure [TSl 
the lifetime could decrease by a factor of 1.2. Likewise, 
the assumed embedded phase lifetime is based on the ra- 
tio of Class and Class I protostars to Class II sources; 
if a large fraction of the Class I source from Evans et al. 
(2008, in prep) are n ot embedded proto stars but edge- 
on TTuari disks (e.g. ICrapsi et al.ll2008l ). the embedded 
phase lifetime could be as low as 2.8 x 10 5 yr, decreasing 
the prestellar core lifetime by a factor of two. 

Finally, the assumption of a simple continuous flow of 
star formation for the last 2 Myr may, of course, be in- 
correct, but without a detailed star formation rate model 
(which might be equally incorrect) it is the best we can 
do. We hope to mitigate errors from variations in the star 
formation rate by averaging over three clouds. Given the 
range of possible errors, we do not apply correction fac- 
tors to the numbers in Tabled but consider the absolute 
uncertainty in the prestellar core lifetime to be a factor 
of two in either direction. 

7. CONCLUSIONS 



cloud density (m easured within the Ay = 2 contour; 



Utilizing large-scale 1.1 mm surveys (|Enoch et al.l 
1200611 Young et al.ll2006UEnoch et afl2007ft together with 
Spitze r IRAC and MIPS maps from the c2d Legacy pro- 
gram (|Evans et al.ll2003h . we have carried out an unbi- 
ased census of prestellar and protostellar cores in the 
Perseus, Serpens, and Ophiuchus molecular clouds. We 
identify a total of 108 starless and 92 protostellar cores 
in the three cloud sample. Based on a comparison of 
1.1 mm derived masses to vi rial masses derived from an 
NH 3 survey of Perseus cores flRosolowskv et al1l2008l ). we 
conclude that the majority of our starless cores are likely 
to be gravitationally bound, and thus prestellar. 

The spatial distributions of both starless and proto- 
stellar cores are similar in these three molecular with 
varying global properties. In all three clouds both star- 
less and protostellar cores are found only at relatively 
high cloud column densities: 75% of cores are associated 
with Ay > 6.5 — 9.5 mag in Perseus, Ay > 6 — 10 mag 
in Serpens, and Ay > 19.5 — 25.5 mag in Ophiuchus. 
Spatial clustering of starless cores is similar in nature to 
protostellar cores but lower in amplitude, based on the 
two-point spatial correlation function and peak surface 
density of cores. 

Cloud environment does appear to have some effect on 
the physical properties of starless cores, however, and 
how they differ from cores that have already formed pro- 
tostars. Starless cores in Perseus are larger and have 
lower mean densities than protostellar cores; we sug- 
gest a simple scenario by which protostellar cores might 
have evolved from starless cores in that cloud, becoming 
smaller and denser at a fixed mass. In Serpens, it ap- 
pears that future star formation will occur in lower mass 
cores than those that are currently forming protostars. 
Meanwhile, in Ophiuchus we see essentially no difference 
between cores that have formed stars and those that have 
not. Of the three clouds, Serpens has the highest mean 



lEnodi et al.ll2007fl and the highest turbulent Mach num- 
ber ([Enoch et alTl2007l ; J. Pineda, personal communica- 
tion), which may be related to its low fraction of starless 
cores (A^sL/^emb = 0.4 compared to 1.0 in the other 
clouds) and the fact that the higher mass cores have al- 
ready formed protostars 

The combined prestellar CMD, which includes 108 
prestellar cores from three clouds, has a slope above our 
completeness limit (0.8M Q ) of a = —2.3 ± 0.4. This re- 
sult is consistent with recent measurements o f the stellar 
initial mass functio n (a = —2.3 to —2.8; e.g. iReid et all 
l2002t lKroupall2002| ). providing further evidence that the 
final masses of stars are directly linked to the core forma- 
tion process. We place a lower limit on the core collapse 
efficiency (the percentage of initial core mass that ends 
up in the final star) of 25%. A more secure link between 
the CMD and IMF requires measurement of the CMD 
down to masses less than 0.2M Q in samples of cores that 
can be demonstrated to be likely prestellar. 

In all three clouds the lifetime of dense prestellar cores 
is similar to the lifetime of embedded protostars, or 
2 — 5x 10 5 yr. The three-cloud average is 4.5±0.8x 10 5 yr 
(with an absolute uncertainty of a factor of two), argu- 
ing strongly for dynamic core evolution on a few free- 
fall timescales. Such a short prestellar core lifetime is 
inconsistent with highly magnetically sub-critical cores, 
in which case evolution should occ ur over a n am bipola r 
diffusion timescale (£ad ~ 10 7 yr; Na kanol e.g.. fl998). 
Our results suggest, rather, a dynamic core evolutionary 
scenario, as might be a ppropriate if turbulence d omi- 
nates the cloud physics (|Mac Low fc Kle ssen 20 0-1), o r 
for near-critical magnetic models (|Ciolek fc Basull2001h . 



The observed prestellar core lifetime decreases with in- 
creasing mean density, at a rate faster than that expected 
from the dependence of the free-fall timescale on density. 

Although this measurement of the prestellar core life- 
time supports a dynamic paradigm over a quasi-static 
one, the distinction is not clear-cut. Our observations 
could still be consistent with a quasi-static picture if we 
are only observing the densest stages (n>2x 10 4 cm~ 3 ) 
of a longer-scale core evolution. Furthermore, the fact 
that we observe extinction thresholds for finding dense 
cores at Ay > 6 mag may be a hint that magnetic fields 
become important in the low column density regions 
of molecular clouds, inhibiting the formation of high- 
density prestellar cores (and thus star formation) . Better 
measurements of magnetic field strengths and more de- 
tailed comparisons between observations and models will 
be help to resolve these ambiguities. 
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Starless cores in Perseus, Serpens, and Ophiuchus 
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NOTE. — Bolocam identification is from Papers I-III. The total mass is calculated from the total flux at 1.1 mm, assuming Tjj — 10 K. 
Deconvolved angular FWHM sizes (#d ec ) an d axis ratios arc derived from an elliptical gaussian fit, deconvolved by the 31" beam. Mean densities 
arc calculated in a fixed linear aperture of diameter 10 4 AU for each source. Note that the Ophiuchus identifications arc different than Tables 1 
and 2 of Paper II for Oph-Bolo 35 and above due to an adjustment of the source list, resulting in 43 rather than 44 sources (IDs arc shifted by -1 
for Oph-Bolo 35 and up). 
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TABLE 5 

Protostellar cores in Perseus, Serpens, and Ophiuchus 
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777 
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Serpens 


Scr-Bolo 2 


18 


28 


44.0 
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0.69 


46.1 


1.2 


1.0x10 s 


Ser-Bolo 3 


18 


28 


45.8 


+00 51 32 


227 


0.52 


35.3 


1.3 


1.1x10 s 


Ser-Bolo 7 


18 


28 


53.0 


+00 19 04 


119 


0.22 


39.0 


1.2 


0.6 x10 s 


Ser-Bolo 8 


18 


28 


55.2 


+00 29 28 


617 


3.55 


59.2 


1.2 


3.2 x10 s 


Ser-Bolo 9 


18 


28 


55.9 


+00 48 30 


141 


0.23 


31.9 


1.6 


0.7x10 s 


Ser-Bolo 13 


18 


29 


00.2 


+00 30 20 


239 


0.48 


44.7 
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1.2x10 s 
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29 


07.0 


+00 30 42 


1016 


2.92 


38.3 


1.1 


4.8x10 s 
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82 


0.54 
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0.4x10 s 
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18 


29 
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337 


2.12 


82.3 


1.0 


1.7x10 s 


Ser-Bolo 22 


18 


29 
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+01 16 51 
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5.10 


39.9 


1.2 


8.1x10 s 
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TABLE 5 









Protostellar 


cores in Perseus, Serpens, 


and Ophiuchus 








Ser-Bolo 23 


18 


29 50.2 


+01 15 


25 


3010 


7.98 


33.2 


1 


.0 


13.2x10 s 


Ser-Bolo 24 


18 


29 53.8 


+00 36 


10 


202 


0.54 


38.0 


1 


.2 


1.0 xlO 5 


Ser-Bolo 25 


18 


29 57.4 


+01 13 


15 


1981 


4.31 


34.1 


1 


.2 


9.5xl0 5 


Ser-Bolo 26 


18 


29 59.2 


+01 14 


07 


1331 


2.56 


41.7 


1 


.2 


6.6xl0 5 


Ser-Bolo 27 


18 


30 00.3 


+01 10 


38 


425 


1.60 


45.5 


1 


.5 


2.2x10 s 


Ser-Bolo 28 


18 


30 00.7 


+01 12 


57 


1266 


2.56 


41.7 


1 


.2 


6.6x10 s 


Ser-Bolo 29 


18 


30 01.0 


+01 


11 


19 


979 


3.30 


40.9 


1 


.3 


4.9 x10 s 


Ser-Bolo 32 


18 


30 05.7 


+00 39 


32 


105 


0.45 


55.8 


1 


.3 


0.5x10 s 


Ser-Bolo 33 


18 


30 06.4 


+00 42 


37 


126 


0.17 


16.7 


1 


.1 


0.4x10 s 


Ophiuchus 


Oph-Bolo 


4 


16 


26 


09.9 


-24 


20 


29 


397 


0.15 


43.0 


1 


.1 


1.3x10 s 


Oph-Bolo 


5 


16 


26 


20.7 


-24 


22 


17 


365 


0.50 


58.5 


1 


.1 


1.2x10 s 


Oph-Bolo 


7 


16 


26 


24.7 


-24 


21 


08 


407 


0.34 


51.2 


1 


.2 


1.2x10 s 


Oph-Bolo 


8 


16 


26 


27.2 


-24 


22 


27 


1382 


1.23 


48.9 


1 


.1 


4.4x10 s 


Oph-Bolo 


9 


16 


26 


27.6 


-24 


23 


37 


2696 


2.02 


38.5 


1 


.2 


8.8 x10 s 


Oph-Bolo 


11 


16 


26 


32.6 


-24 


24 


45 


1249 


0.54 


39.8 


1 


.1 


4.0 x10 s 


Oph-Bolo 


12 


16 


27 


00.7 


-24 


34 


17 


540 


1.74 


109.8 


1 


.2 


1.9 x10 s 


Oph-Bolo 


14 


16 


27 


07.9 


-24 


36 


54 


261 


0.37 


76.9 


1 


.6 


0.6x10 s 


Oph-Bolo 


17 


16 


27 


22.3 


-24 


27 


36 


404 


0.32 


52.2 


1 


.2 


1.2x10 s 


Oph-Bolo 


18 


16 


27 


25.2 


-24 


40 


29 


172 


1.17 


114.6 


1 


.6 


1.3x10 s 


Oph-Bolo 


19 


16 


27 


27.0 


-24 


26 


57 


615 


0.25 


45.5 


1 


.1 


2.3 x10 s 


Oph-Bolo 


20 


16 


27 


29.1 


-24 


27 


11 


691 


0.27 


42.9 


1 


.1 


2.3 x10 s 


Oph-Bolo 


26 


16 


28 


21.0 


-24 


36 


00 


225 


0.51 


114.8 


1 


.2 


0.6x10 s 


Oph-Bolo 


30 


16 


31 


37.2 


-24 


01 


52 


289 


0.14 


39.0 


1 


.2 


0.5x10 s 


Oph-Bolo 


33 


16 


31 


52.6 


-24 


58 


01 


204 


0.18 


49.3 


1 


.2 


0.5x10 s 


Oph-Bolo 


35 


16 


32 


00.6 


-24 


56 


11 


189 


0.14 


63.3 


1 


.2 


0.4 x10 s 


Oph-Bolo 


36 


16 


32 


24.3 


-24 


28 


19 


4624 


1.16 


17.6 


1 


.3 


4.7x10 s 



Note. — Table columns arc as for Tablc[4] Note that the Ophiuchus identifications arc different than Tables 1 and 2 of Paper II for Oph-Bolo 35 
and above due to an adjustment of the source list, resulting in 43 rather than 44 sources. The peak flux, total mass, sizes, etc. for Oph-Bolo 35 
and Oph-Bolo 36 arc also different than in Paper II. 



